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Abstract. - We report on the evolution of the Raman-active A g phonon modes in the charge- 
ordered manganite Pro.63Cao.37Mn03 as a function of temperature from 300K to 25K. Our 
studies reveal that the linewidths of the A 9 (2) and A 9 (4) phonons increase as temperature 
decreases. This anomalous temperature dependence of phonon lineshapes, seen for the first 
time in charge ordered manganites, can be quantitatively understood in terms of a strong 
spin-phonon coupling involving t^g spins and A g phonons. 



The recent upsurge of interest in doped manganites exhibiting colossal magnetoresistance 
has fuelled a lot of interesting questions regarding the role of spin, orbital and lattice degrees of 
freedom and the interplay amongst them Q. These systems of the form Ai_ x A' x Mii03, where 
A is a trivalent rare earth ion (e.g., La, Pr, Nd) and A' is a divalent ion (e.g., Ca, Sr, Pb) show 
a rich phase diagram depending on the tolerance factor and the amount of doping x which in 
turn controls the ratio of Mn 3+ to Mn 4+ . The essential ingredients for understanding their 
electronic properties are the double-exchange (DE) mechanism J| and the polaron formation 
due to the Jahn-Teller(JT) effect |]. In the case of systems where the weighted average 
A-site cation radius is small, the system also shows charge ordering (CO), i.e., real space 
ordering of Mn 3+ and Mn 4+ ions at low temperatures [0. The CO state becomes stable when 
the repulsive Coulomb interaction between the carriers or the cooperative JT effect dominates 
over the kinetic energy of carriers. In such situations there is a strong competition between the 
DE interaction which favours ferromagetism (FM) , and CO which favours antiferromagnetism 
(AF) via the coupling between the background t^g (S = 3/2) spins. 
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The system under study, Pri-^Ca^MnOs, shows a rich phase diagram as a function of 
doping and temperature and has been studied using a variety of probes like resistivity j5|, 
magnetization and neutron diffraction ||. In the temperature range T > Tco (~ 240 K) the 
system is a paramagnetic insulator. In dc magnetic susceptibility, a large peak is observed at 
Tco followed by a relatively small peak at Tat. The peak at Tco is attributed to ferromagnetic 
correlations [Q. The system further undergoes a transition at Tat ~ 17 5K to the charge- 
exchange (CE) type antiferromagnetic insulator (AFI). The CE type structure for x < 0.5 
is different from what is present in systems with x = 0.5. In this so called "pseudo-CE" 
structure the zig-zag FM chains in the ab plane are FM aligned along the c-axis [[?J. This 
is unlike the CE structure where the layers in the ab plane are aligned antiferromagnetically 
along c. Further lowering of temperature below 50 K leads to another transition which can 
be understood either as a canted antiferromagnetic state || or as a mixture of ferromagnetic 
domains (or clusters) in an antiferromagnetic background. 

In this letter we report our Raman studies on well characterized, high-quality single crystals 
of Pro.63Cao.37Mn03 ||] as a function of temperature from 300K to 25K. There are many Ra- 
man studies on manganites in recent times covering both the vibrational [p-13 and electronic 
Raman scattering |14| across the metal-insulator transition. In a recent work on polycrys- 
tallinc Pro.65Cao.35Mn03 temperature dependence of the modes at 475 and 610cm _1 
was studied and the splitting of these modes below Tco was attributed to the folding of the 
Brillouin zone in the CO state. In another Raman study Jl6|], the effect of A-site substitution 
on the A s (2) mode was discussed. However, the linewidths have not been discussed, and as 
we show in this letter, they reveal a very unusual behaviour. The aim of the present study of 
P r o.63Cao.37Mn03 is to explore the coupling between phonons, electrons and spins by examin- 
ing the temperature dependence of Raman linewidths across the CO and magnetic transitions. 
Our studies do reveal a strong anomalous temperature dependence of the linewidths of the two 
Raman active A g modes, which we have quantitatively analysed in terms of strong spin-lattice 
coupling in these systems. 

The single crystal used for the experiment was prepared by the float zone technique Q. 
A (100) crystalline face of dimensions 5mm diameter and 2mm thickness was polished using 
diamond paste and mounted on the cold finger of the closed-cycle helium refrigerator (RMC 
model 22C CRYODYNE) using thermally cycled GE (M/s. General Electric, USA) varnish. 
Raman spectra were recorded in the spectral range 150 cm -1 to 750 cm -1 in back-scattering 
geometry using a DILOR XY spectrometer equipped with a liquid-nitrogen-cooled charge- 
coupled device detector. An argon ion laser line at 514. 5nm was used as the excitation 
wavelength at laser power of ~ 25mW focussed to a spot of 40/im on the sample. The 
temperatures quoted are those measured on the cold finger using a Pt-sensor coupled to a 
home-made temperature controller. The temperatures are within an accuracy of ± IK and 
do not take into account laser heating. 

For an ideal perovskite ABO3 with cubic structure (space group Ol, Pm3m) the irre- 
ducible representations at the zone center are given by 4Fi„ + F2 U - The only allowed optical 
modes are the infrared modes and since all lattice sites have inversion symmetry, first-order 
Raman scattering is forbidden. The small size of the Pr cation in comparison to the oxygen 
causes a large tilting of the MnOg octahedra leading to a GdFeOa structure. The space group 
of Pr .63Cao.37Mn03 at room temperature is Pbnm (D^); there are four formula units per 
unit cell and there are 24 Raman active modes || . Iliev et al. ^ have carried out systematic 
Raman measurements at room temperature and lattice dynamical calculations on undoped, 
stoichiometric, orthorhombic YMn03 and LaMnOs with Pbnm structure and identified the 
symmetry of the different modes. The Raman modes are weak even in stoichiometric LaMnOs; 
the observed modes at 140, 198, 257, 284 and 493cm~ 1 were assigned A g symmetry, modes 



R. Gupta el al: Anomalous Raman Linewidths in CO Manganites 



3 




Fig. 1 Fig. 2 

Fig. 1 - Raman spectra as a function of temperature, (a) Micro-Raman spectrum taken at 298 K 
showing the modes at ~ 260, 450, and 650 cm -1 , (b) Raman spectrum of A 9 (2) and A 9 (4) modes at 
three different temperatures. The thick solid smooth curves show the sum of the two fitted Lorentzian 
lineshapes (thin solid smooth curve). Raw data has been denoised using a wavelet filtering routine 
(Daubechies wavelet filter of order 4 at level 4). 

Fig. 2 - Variation of the Raman frequency (lu) and the full width at half-maximum (F) as a function 
of temperature. Filled squares are the experimental values and thick solid curves are the theoretical 
results in our model calculation. Panel (a) is for the A 9 (2) mode and panel (b) is for the A s (4) mode. 



at 109, 170, 308, 481 and 611cm -1 as B 2g and modes at 184 and 320cm -1 as B lg or B 3g . Re- 
cent experiments on LaMnOa have revealed the presence of high-energy modes (~ 1000cm -1 ) 
which have been attributed to the orbital excitations ]l7j| . However, doped manganites be- 
ing pseudo-cubic, Raman lines are extremely weak and broad |^13]| . The number of modes 
observed in doped manganites are very few; e.g., in Lao.7Cao.3Mn03 there are five modes 
observed at low temperature at 70, 133, 235, 435, and 670cm -1 O]. 

In the micro-Raman setup using 6328 Aradiation from He-Ne laser, we observe four bands 
centered around ~ 260, 290, 450 and 650cm -1 . However, in macro-geometry with the crystal 
in the cryostat, only two modes centered around 258 and 288cm -1 at 300K were observed 
with sufficient signal to noise ratio (Fig. 1). The other bands near 450 and 650cm -1 are 
weaker and therefore could not be followed as a function of temperature. 

We will now discuss the origin of the observed modes. Since the crystal structure for 
Pr 1 _ 2 ,Ca 2 .Mn03 is same as that of LaMn0 3 , it is plausible to assign to these modes the 
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symmetry used by Iliev et al. j|] for LaMnC>3. The low-frequency mode at 258cm" 1 has been 
attributed to the A ff (2) mode which is the in-phase rotation of the oxygen cage about the y-axis 
on adjacent MnOg octahedra. In this mode only the in-plane oxygen atoms, namely 01, are 
involved. The mode at 288cm -1 has been assigned to the A ff (4) mode which is out-of-phase 
rotation of the oxygen cage about the x-axis. Here along with the in-plane oxygens (01), the 
apical oxygens (02) are also involved. The Mn ion being at the centre of inversion symmetry 
remains stationary in both the modes. Figs. 2a and 2b show the variation of the peak position 
and lincwidths (full width at half-maximum) for the two modes as a function of temperature, 
obtained by nonlinear least-square fitting of the data with a sum of two Lorentzians along 
with a baseline. The mode frequencies increase as the temperature is lowered and a sharp 
change is noticeable at ~ 50 K as the system undergoes a transition to a canted AFI phase Q . 
Most interestingly, the linewidths for both the modes show anomalous behaviour in that the 
widths decrease as the temperature is increased. Another notable observation is the large 
value of linewidth even at low temperatures. 

Next we look at the various possibilities that may lead to the observed temperature de- 
pendence of mode frequencies. Mode softening (~ 3%) as the temperature is raised can arise 
due to the quasi-harmonic contribution arising from the lattice expansion which can, in turn, 
change the force constants. This change in frequency is related to the Griineisen parameter 
and the volume change. Since the fractional volume change across the studied temperature 
regime is not available for Pro.63Cao.37Mn03 we can be guided by the data available in liter- 
ature for Pro.7Cao.3Mn03. The volume change is ~ 0.4% between temperatures of 300 and 
10 K. This change is too small to account for the observed frequency changes of both the 
A g modes. Another contribution can arise due to the intrinsic anharmonic (cubic) interac- 
tions. This contribution is also insufficient to explain the large softening, since typically for 
perovskites, e.g., KTa03, the change of peak position is about 2-3 cm" 1 in the temperature 
range 10-500 K @. 

Now we consider the linewidths. Anharmonic interactions will lead to an increase of 
the linewidth with increasing temperature, in contrast to our experimental results. Lee and 
Min jl9| have proposed a model to explain the sound velocity softening across the charge 
ordering transition in Lao.sCao.sMnOa pCfl , Recently Dattagupta and Sood Q have extended 
this model incorporating both the DE interaction and electron-optical phonon coupling. Their 
models also fail to explain the behaviour of linewidth at low temperatures in the magnetically 
ordered phase. We suggest that magnon-phonon interaction leads to the shifts in phonon 
frequencies as well as the observed temperature dependence of lifetimes. The frequency shift 
can be understood as arising due to the modulation of the exchange integral between Mn ions. 
This has been considered earlier for LaMn03+^ 13 to explain the softening of the Mn-O-Mn 
stretching mode. However, this model also does not take into account the dynamics and 
hence cannot explain the temperature dependence of the linewidth(r). A sharp increase in 
the peak positions at low temperatures near the canted AF transition (see Fig. (2)) suggests 
a strong coupling between the Raman modes and spin excitations, which we show to be also 
responsible for the anomalous temperature dependence of the linewidths. 

We briefly outline here the origin of the large spin-phonon coupling and the way it leads to 
anomalous phonon damping. Details will be published elsewhere ]22] ]. When the neighbouring 
MnOg octahedra are tilted, the Mn-O-Mn bond angle %(= [tt — (Oi — Oj)}) between nearest- 
neighbour Mn-0 bonds involving Mn ions i and j changes. Here Oi and Oj are the small 
angles of rotation of the neighbouring octahedra centred at i and j, respectively. Since the 
oxygen mediated overlap t%j between Mn-d orbitals goes as cos Oij, Jij goes as Jfj cos Oij/ 2 ~ 

Jij(l — 2 )• ^his nas a term of the form Qfij which, among others, corresponds to a 
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simultaneous excitation of the rotation mode involving the Mn-0 bond, with a displacement 
adi) and another mode q2j (~ a9j) where a is the Mn-0 bond length. The coupling 
between the two phonons will be strongest if they belong to the same symmetry representation. 
The first of these is the Raman-active phonon under consideration, with a frequency wi = 
258cm -1 (A g (2)) or uo\ — 288cm -1 (A g (A)) and the other is an optical phonon mode with 
frequency u>2- 

For simplicity in calculation, we assume that the magnetic order is of the C-type in which 
ferromagnetic chains are linear in contrast to the zig-zag chains in the CE structure. Both 
these structures are geometrically similar, have the same number of nearest-neighbour FM 
and AF bonds per spin, and have the same magnetic energies per spin p3[ . Though there 
can be quantitative differences in their magnetic excitation spectrum, qualitatively both are 
very similar in that they have mixed FM and AF character along arbitrary directions of the 
Brillouin zone with the same number of nearest-neighbour FM/AF couplings and the same 
overall geometry, namely FM along a line and AF in plane. The latter decides the dispersion 
relation and the phase space for the mixed FM and AF spin excitations. Since we are interested 
in understanding how magnetic fluctuations affect the phonon spectrum, it is therefore a very 
good approximation to work with the C-type structure rather than the CE or pseudo-CE type 
structures. The effective spin-Hamiltonian for this anisotropic spin system is given by 



where AF and F refer to antiferromagnetic and ferromagnetic interactions. Our Hamiltonian 
is most naturally viewed as an effective Hamiltonian with t^g spin degrees of freedom kept 
explicitly, and all other degrees of freedom (such as e g electrons, JT coupling between e g elec- 
trons and lattice modes orthogonal to those considered here, Hund's rule coupling) integrated 
out. (This is possible and sensible because the t 2g spins are the lowest energy degrees of 
freedom near and below Tat, as discussed below.) It has been shown by several authors, e.g., 
Kugel and Khomskii |m| , that in LaMnOs (which is an insulator experimentally showing JT 
distortions and ferromagnetic planes coupled antiferromagnetically) the Mn-O-Mn-mediated 
magnetic interaction is /erromagnetic in plane and antiferromagnetic (usual superexchange) 
perpendicular to it. The latter is indeed due to the same sort of virtual processes involving the 
Mn-O-Mn bond as the AF superexchange. It has been shown that [Q the JT distortions (and 
consequent e g orbital excited states) , and the Hund's rule coupling Jh (which enslaves the e g 
spins to the t2 g spin in the ground state and produces an excited state an energy Jh away, 
with e g , t2 g spins antiparallel) lead to a ferromagnetic Mn-O-Mn-mediated in-plane exchange 
coupling. Thus detailed models involving DE and JT coupling show as briefly mentioned 
above, that both AF and F coupling between t2 g spins necessarily involve Mn-O-Mn overlap. 
Further, if there are no other low-energy degrees of freedom (e.g., no e g fermionic low-energy 
metallic excitations, or orbital fluctuations), then the form of Eq. (1) of the text with the 
discussed dependence on Mn-O-Mn angle, is inevitable. A microscopic theory will relate the 
Jij to basic parameters of the electron Hamiltonian, but Jy is always second order in the 
Mn-O-Mn overlap. The only other possible low-energy degrees of freedom in the insulating 
phases of CMR oxides are orbital fluctuations. Here we contend that these are not relevant 
in the temperature range we are working in because orbital order is already well developed. 
For small fluctuations of the Mn-0 bond angle, from the above discussion, the spin-phonon 
interaction Hamiltonian (H s - P ) becomes 




(1) 
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+ JoAF^i ' Sj ^ ^oF^i ' Sj, (2) 

(ij);i,jexz (ij);j£y 

where 5 is a coupling constant of order unity. The spin degrees of freedom are described in 
terms of Schwinger bosons (SB) p5| which is an exact representation. It is convenient both 
below and above the Neel temperature (TV). Below T^r, the SB Bose condense. The spin-spin 
interaction in Eq. (2) can be written in terms of boson operators a ka in mean-field theory as 

Hsb = Hp + }^k[a krT aho + !]• (3) 

Here Hq is a c-number (with no operators in it unlike the second term which is the Hamiltonian 
of noninteracting SB) to be determined self-consistently from the free energy, and u>k is the 
boson frequency, 

= -y/p - J y oF z y B yl l) + A] 2 - (J*Z F z xz B xzl r) 2 - (4) 

Here, B xz and B y are the Weiss fields (the SB order parameters), and A is the temperature 
dependent stiffness of the SB excitations, to be determined self-consistently from the free 
energy; z xz (z y ) is the xz plane(y-axis) coordination number, and 7 p = l/z p ^2 S e lkp ' S (p = xz 
or y). We choose J*af ~ 3.9meV and J% F ~ 7.8 meV so as to give T^ = 175 K. 

The spin-phonon interaction in Eq. (2), written in terms of the SB operators, leads to 
a coupling between the A g mode being observed, and the second optical phonon mode of 
wavenumber LO2, with the excitation of nearest-neighbour spin fluctuations or boson pairs. 
We calculate the decay rate for the A g phonons of frequency Tlu)\ in perturbation theory, to 
second order in H s - p . The phonon decays due to the second-order process which involves the 
A g phonon decaying into the second phonon and two Schwinger bosons. The imaginary (real) 
part of the self-energy E determines the decay rate (peak shift) of the Raman-active phonon. 
Assuming the phonon modes to be dispcrsionless, the imaginary part of E is given by 

Im E = - — j 2j/(g,fc)(J(wi -Lo 2 -Wfc - u g )[l + g u ( q+k ) + g^i 1 + 9w{k) + 9w(q)], (5) 

q,k 

where /(g, £;), the form factor is given by 

/(<?,*) = (J! F B y ) 2 p 1 Al + (jy F B y )(J^ F B xz )A 1 (p 2 A 2 +p 3 A 3 )/2 

+ (J^ F B xz ) 2 (p 4 Al+p 5 A 2 2 + P6 A 2 A 3 )/4, (6) 

Pi are products of Bogoliubov or coherence factors, g~} k \ = e /3uk — 1, A% = 2(cosq a + cosfc a ), 
A 2 = 2(cosg x + cosgz), and A 3 = 2(cosfca; + cosfc z ). 

We find substantial linewidth at T — 0, which initially increases with temperature because 
of thermal availability of two SB states at energy Ti{w\ — lo 2 ). At high temperatures, of order 
Tjv or higher, the spin fluctuation mode softens, and the density of pair excitations at a 
fixed energy difference K(u)i — uj 2 ) decreases as temperature increases. The Weiss fields also 
decrease with increasing temperature. Thus the A g (2) phonon decay rate falls. This decrease 
has been estimated in the SB formalism both below and above T^r (the Bose condensation 
temperature in that formalism) , or in a classical theory for spin dynamics of a Heisenberg spin 
system specially useful above Tjv; the results above T^v arc roughly the same. In the latter 
theory, one essentially obtains the nearest-neighbour spin-fluctuation frequency spectrum via 
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classical equations of motion of coupled spin vectors. We show the Schwinger boson results 
in Fig. 2 for both linewidths and peak positions. Guided by the lattice dynamic studies of 
Iliev et al. |j] on LaMnO-3, we have taken the A g phonon at 190cm -1 as u>2- The qualitative 
behaviour of T and peak postions will not change with another choice of u>2- As can be seen, 
the agreement between the experiment and theory is rather good above 70K. Our approach 
does not take into account the observed canted AFI state below T < 50K, which may be 
responsible for the discrepancy between the theory and experiments at low temperatures. 
The canting of spins leads to FM correlations in an otherwise perfectly ordered AF state 
which, in turn, modifies the spin excitation spectrum. 

In conclusion, our Raman measurements on Pri_j.Caj.MnO3 (x = 0.37) reveal a strong 
coupling between the spin and lattice degrees of freedom in CO manganites. The large vari- 
ation of the mode frequency and anomalous behaviour of the linewidths as a function of 
temperature are understood in terms of spin-phonon interaction. Thus, we believe, our study 
reveals the striking consequences of the strong interplay between the spin and lattice degrees 
of freedom in CO manganites. 

RG and GVP thank the Council of Scientific and Industrial Research (India) and AKS 
thanks the Department of Science and Technology (India) for financial support. 
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